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Abstract 



We analyze here a model for an adsorbate system composed of many layers by extending 
a theoretical approach used to describe pattern formation on a monolayer of adsorbates 
with lateral interactions. The approach shows, in addition to a first order phase transition 
in the first layer, a transition in the second layer together with evidence of a "cascade" of 
transitions if more layers are included. The transition profiles, showing a staircase structure, 
corroborate this picture. The adsorption isotherm that came out of this approach is in 
^v^j ' qualitative agreement with numerical and experimental results. 
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1 Introduction 



The problem of phase transitions and ordering phenomena in adsorbed films on crystal surfaces 
has attracted considerable interest for many years (see for a recent review). Formation of 
ordered structures of particles adsorbed on surfaces due to mutual (lateral) interactions has been 
observed in many experiments, originating different types of adsorption isotherm g, [|. Strong 
attractive (lateral) interactions lead to phase separation with different coverages, characteristic 
of first order phase transitions j|, |3| . 

It is well known experimentally that adsorption on surfaces is not only restricted to the 
formation of monolayers, but rather a second layer can condensate on the first one, a third on the 
second, and so on Numerical simulations also show similar results ||. The old BET theory 
|| considered the possibility of multilayer adsorption but, as it does not include interactions 
among particles in the same layer, does not yield either a phase transition or ordered structures. 
So far we are not aware of theoretical models to describe this situation. Other theoretical forms 
for isotherm intended to describe such a multilayer adsorption phenomenon have been discussed 
in the literature p]]. However, such models do not seem to be able to describe phase transitions 
or pattern formation. 

The aim of the present work is to introduce a simple model that can describe multilayer 
adsorption qualitatively and the formation of ordered structures. To reach this goal we exploit a 
model introduced in Ref.(2| to describe a system of a monolayer of adsorbed particles with lateral 
attractive interactions leading to the formation of interfaces separating phases with different local 
coverage. The new aspect is that we extend such a model and go beyond the one layer system 
considering a multilayer system (however, focusing on the case of only two layers) including 
also an exchange interaction between different layers. It is worth remarking here that such a 
"macroscopic" model can be derived from a microscopic one following a procedure similar to 
that used in Q. 

In what follows we present the model, show the form of the isotherm, and discuss briefly 
some aspects of the critical point and the form of the density profiles of the possible fronts. 
Finally we draw some conclusions. 

2 Adsorptive multilayer system 

We adopt here a phenomenological point of view. However, a formal derivation of the macro- 
scopic master equation leading us to the phenomenological equations we present here, is included 
in the appendix, where we start from a microscopic description closely following the procedure 
of Ref.g). 

We adopt a continuous description for the surface, and characterize the adsorptive species 
through evolution equations for q(x, t), the local coverage in the i-th layer (see the appendix for 
the precise definition of these local coverage). The adsorptive contribution is characterized by 
k a (a constant), and the adsorption is only possible in the (1 — q) free sites of the corresponding 
layer. Hence, the adsorption rate is k a p(l — q), where p is the partial pressure of the gaseous 
phase. Additionally, for adsorption on the i-th. layer, it is required to have free sites on the 
i + 1-th layer. 

The desorption process has a rate kd, that includes kd t o, the desorption for noninter acting 
particles, and a correction due to the lateral interactions. The strong local bond induced by 
the interaction Ui(x) (U2(x) for the second layer, etc), corrects the desorption rate as: ka = 
kd t o exp[Ui(x) /kT], where k is the Boltzmann constant and T the temperature. According to the 
form we use to introduce such an interaction, we are assuming that it is a substratum mediated 
interaction. A desorption process in the first layer can happen only if the corresponding site on 
the second layer is free. 
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The potential Ui(x) produces a force Fi = — 9U g^ , that affects the adsorbed particles in- 
ducing a velocity Vi = bFi, where b is the mobility (given by Einstein's relation b = with D 
the diffusion coefficient in the i-th layer) . The associated particle flux is ji ~ c% v $ . Because the 
flux is only possible at the (1 — Cj) free sites, its final form is 

D dUi{x) 



KT l/ dx ' 

while the diffusion flux is given by 

jdif = _ D ^_ jdif _ D dc 2^ 
1 dx ' 2 dx' 

for the first and the second layer respectively. 

The possible transference interaction between different layers is characterized by a coupling 
parameter kj*. It involves free sites in the final layer and occupied ones in the original layer. 
Hence, the form of this contribution is fcrCj(l — q+i) for transfer from a given layer to the upper 
one, and similarly in other cases. In Fig. 1 we sketch the different processes entering into our 
model. 

Considering all the above, the evolution equations for the different coverage can be written 

as 

d 

ci(x,t) = k a p(l - ci)(l - c 2 ) - k dfi (l - c 2 )c 1 ex.p[Ui(x)/kT] + a(l - c 2 )k T 



dt 



d 

c 2 (x,t) = k' a pci(l - c 2 ) - k' d0 C2exp[U 2 (x)/kT} + k' T c 2 (l ~ ci) 



dt 



d r D dU 2 (x) , . ^dc 2 , , \ 



According to [||, we can employ the following simple approach. In thermal equilibrium the 
coverage c%(x) of the i-th layer is stationary, that is ^ = 0. The flux induced by the gradient 
of the potential is balanced by the diffusion one yielding 

k a p(l- Cl )(l-c 2 ) k d>0 (l - c 2 ) Cl e u ^ kT + d(l - c 2 )k T = (3) 
k' a p Cl (l-c 2 ) - k' dfi c 2 e u ^' kT + k' T c 2 {l- Cl ) =0. (4) 

If we consider the case of constant Ui{x) as well as no transfer between the layers, we have 

k a p{l- Cl )(l-c 2 ) - k d>0 {l - c 2 ) Cl e u ^ kT = (5) 
k' aPCl (l-c 2 ) - k' d>0 c 2 e u ^ kT = 0. (6) 

We can get c\ from the first expression obtaining a result that agrees with the well known form 
of the Langmuir isotherm for adsorption in a monolayer |§] 

ci = • (7) 

1 + (& Ul/kT 

We can also we obtain c 2 from the previous system and consider the total coverage, that can be 
eventually compared with the results of the usual BET isotherm || or similar ones 0. However, 
our main interest is to obtain relations and conditions yielding the phase transition, as well as 
the form of evolution fronts. 
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For the functional form of Ui(x) we assume an attractive (and as indicated earlier, substratum 
mediated) potential among particles separated a distance r, that we denote by Uj(r). The 
potential acting on a particle located at r in the first layer is 

U\(r) = — / ui(r — r')c(r') dr' '. (8) 



The integration domain is the whole surface. As a first step in our analysis we assume that we 
have the same potential u(r — r') among particles in all layers. The function u(r) depends on 
the nature of the system. If the interaction radius is small compared with the diffusion length, 
and the coverage is not much affected by variations in this radius, we can approximate 

/dc 1 d 2 c (r^ 
u(r - r')h(r) + (r - r')—±- + -lr - r') 2 ' H ]dr'. (9) 

or 2 or A 

The spatial derivatives are evaluated at r. Hence we obtain 

f d 2 c\ 
u i( r ) = - u(r - r')ci(r')dr' ~ -UqC\ - Xl-^i ( 10 ) 

where the coefficients are given by 

no = J u(r)dr, (11) 

and 

Xi = 1/2 / r 2 u(r)dr, (12) 



and from symmetry considerations we have / ru(r)dr = 0. 

In a similar way (and remembering that we assume the lateral interaction is substratum 
mediated) for the second layer we get 

U 2 (r) = ci[- / u(r - r')c 2 {r')dr'}. (13) 



After replacing C2(r') by its approximation (that is a Taylor expansion around r as in Eq. @), 
we find 

f d 2 C2 
U 2 (r) = -ci / u(r - r')c 2 (r')dr' ~ -ci[u c 2 + X2^4], (14) 



dr 2 

where no and \2 have the same form as before. 

Replacing U\{r) y U 2 (r) into the equations of motion Eqs.(l,2), we get 

9 Cl (r,t) = k aP (l - Cl )(l - c 2 ) - k dt0 {l - c 2 ) Cl e- U0C1 / kT + ci(l - c 2 )k T 



Ot 



+ frW-wt 1 -*™^ (15) 



9 P2(r,t) = k' a p Cl {\ - c 2 ) k' dfi c 2 e~ c ^l kT + k' T c 2 (l - d) 



+ ^1-^(1-^1. (16) 



In both equations we have neglected the contribution coming from the terms Xi9 2 Ci, as we follow 
Ref. H, and assume they are only small contributions. 



4 



To simplify the notation we assume D\ = Z? 2 = D, and we scale the variables as follows: 
£ = x/Ljif, where the diffusion length is L^/ = (D/k^o) 1 ^ 2 ; r = t/td, with = l/fc^o! 
e = uo/kT. Finally, a = k a p/kdfl characterizes the coverage in equilibrium when U(x) becomes 
zero. Also, K = kx/kdfi is the exchange parameter. 

With the indicated scaling the evolution equations for Cj(£, r) reduce to 

— = ai(l - ci)(l - c 2 ) - (1 - c 2 )ci exp[-eci] 

- Cl )^] + ^ + d(l - c 2 )K (17) 



9 ° 2 n \ r 1 

— = a 2 ci(l - c 2 ) - c 2 exp[-ecic 2 J 



- c 2 )c 2 ^] + ^* + c 2 (l - d)X. (18) 



As in 0, we may consider for the stationary regimen or thermodynamic equilibrium (and 
initially assuming K = 0) that 

«i(l — ci) — ci exp[— eci] = (19) 
a 2 ci(l - c 2 ) - c 2 exp[-ecic 2 ] = 0. (20) 

The dependence of the coverages q on e, «i and a 2 is apparent. The result is analogous to 
that shown in Fig. 2 of Ref.Q. It shows that in the plane (e,cti) there are three regions. Two of 
them correspond to real solutions, with homogeneous coverage, while the intermediate region is 
bistable and presents a coexistence phase with intermediate coverages. Those regions coalesce 
in a critical point indicating a second order phase transition. 

In equilibrium and from the diffusive term we get equations for c\ and c 2 : 

1 — eci(l — c\) = 0, 
1 - eci(l - c 2 )c 2 = 0. 

The last expressions indicate that in order to have a high coverage in the first layer we need 
e > 4, while for the second layer we find a relation with the coverage of the first one through 
e > 4/ci. This indicates that c 2 becomes relevant only after some threshold value of the coverage 
has been reached in the first layer. 

The kind analysis we have just outlined, has been strongly criticized in on the basis that 
the Maxwell construction in [0] is done assuming that coverages near the interface correspond 



to the spinodal points. This was discussed and corrected in [1C] using the true equilibrium 
condition, that is that the chemical potentials of both phases should be equal at the transition 
point. In the next Section we use such an approach in order to determine the correct phase 
diagram. 

3 Local Chemical Potential 

As indicated in the last part of the previous section and in order to analyze the phase transition 
we follow |L"o| l defining a local chemical potential for each layer according to 

m = - [u(r- r') ci (/) dr' + kT ln[ Cl ^\ - 1 , (21) 
J 1-cAr) 



5 



p' = -ci / u(r - r')c 2 {r')dr' + kT Inf ° 2 , J - (22) 
7 I - coir) 



C 2(r) 
c 2 (r) J 

Using these definitions, the equations of motion can be written as 



where 



and 



^ = M l-c 1 )[l-^]+^[^c 1 (l-c 1 )^], (23) 
~i = k aP ci(l ~ c 2 )[l - e^] + ^c 2 (l - <*)-£], (24) 

Po = *Tln[^], (25) 

^ = fcr ln[^ Cl ]. (26) 

In equilibrium these chemical potentials are constant on all surface regions corresponding to 
each layer and equal to <po, for the first one and cp' Q , for the second. 

For a uniform distribution, according to Eqs.(^ll) and (|22|), the chemical potentials are given 

by 

tp(tn) = -u ci + kT ln[-^] , (27) 

1 - Cl 

and 

v '(c 2 ) = - UoCl c 2 + kT ln[-^-]. (28) 

1 — C2 

In figure 2 we depict the phase diagram that clearly shows three regions: a low and a high 
coverage region, and a central region corresponding to coexistence. This last region is bistable 
in both layers and fulfills the conditions = 0; = 0; that explicitly indicate 

dip kT 

7^ = -^ + —p. r = 0, 29 

and 

dip' kT 

-f- = -u QCl + — = 0. 30 

dc 2 c 2 (l - c 2 ) 

In order that the interface separating both phases be stationary, the Maxwell condition 
should be satisfied in both cases. Such a condition can be derived for p = po and p' = p' , 
together with the approximation indicated in Eq. (|io|) . This yields the following equations 

ci <9 2 ci 

- uqC\ + kT ln(- )- Xl —— = ip 0: (31) 

1 — ci or z 

and 

c 2 d^c 2 

- u c 2 + kT \n(-^——) - X 2-^2 = ¥o- (32) 

^From these expressions we can obtain the Maxwell condition corresponding to the coexistence 
of both phases. For the first layer we have 

fCi2 

(<p(ci) - (po)dci = 0, (33) 

'en 

where ip{c\) is given by Eq. (p7|), and en, ci 2 are the coverage values in the equilibrium phases 
on the first layer. For the second layer we have 

r c 22 

(p'{c 2 ) - <p' ) dc 2 = 0, (34) 

'C 2 1 

where <p(c 2 ) is given by Eq. (p8|), and c 2 i, c 22 are the coverage values in the equilibrium phases. 
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4 Stationary coexisting states 
4.1 First layer 

As seen above, the coexistence line between the dense and diluted phases on the first layer is 
given by the Maxwell condition in Eq. (|33"|), where 



(p(ci) = -u ci + kT ln( ^ ), (35) 

1 C\ 

= kT In — ^ (36) 
kd,o 



After integrating we find 



2 

^° Cl1 + kT ln(l - en) + kT cu ln( 11 ) - yocii = 

2 1 — Cu 

= ^^ + A;rin(l- Cl2 ) + A ; r Cl2 ln( T ^^)-^oCi 2 . (37) 

2 1 - C12 

Here, cu y ci 2 are the equilibrium coverages of both phases on the first layer. They satisfy the 
equations 

- u c n + kT ln( ° n ) = v? , (38) 

1 - cu 

- noci2 + feT ln( ° 12 ) = (fo, (39) 

1 - C\2 



with given by 



cp = kT In(^E). (40) 



After some algebra, where the condition cu + ci 2 = 1 arises, we get 



- u c n + kT ln( Cn ) = --tt - (41) 
1 - cu 2 

When comparing Eqs. ( |38| ) and (|4"l"D, we see that 

¥>o = ~^ n o- (42) 
Replacing ifQ by Eq. (|4(]) we finally obtain the coexistence condition on the first layer 

k a p __^L tAO , 
- — = e 2kT . (43) 

kd,o 

4.2 Second layer 

In a similarly way for the second layer, we consider Eq. (|34|), where 

¥>' = -« cic 2 + kT In ( 3-^-) , (44) 

and 

( h'm \ 

(45) 
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After integrating we get 

- U ° ClC21 + kT ln(l - c 2 i) + kTc 2 i ln( ° 21 ) - ip' c 2 i = 
I 1 — C21 

= -^^ + A;rin(l-c 22 ) + A ; Tc 22 ln(-^-)-^c 22 . (46) 
z 1 — C22 

Here c 2 i and c 22 are the equilibrium coverages in both phases of the second layer. They fulfill 
the equations 

- n c lC21 + kT ln(-^— ) = v ' , (47) 

t - c 2 i 

and 

- u C!C22 + kT ln(- ) = ip' , (48) 

1 - c 22 

^ = kT ln(^ Cl ). (49) 

As in the previous subsection, after some algebra where the condition c 2 i + c 22 = 1 arises, 
we get 

— UqC\C2\ + kT ln(- — — — ) = —-UqCi. (50) 
1 - c 2 i 2 



and ip' is given by 



Comparing Eqs. ( [471) and (|50|), we see that 



Replacing by Eq. ((f49|) we get 

For the first layer we denote 
while for the second 



.1 

(f = --u ci. (51) 



k' a p 1 _H2£1 ,_ nN 

-^il = — e 2kT . (52) 

Kd,0 Cl 



5 Adsorption isotherm and Critical behaviour 



The adsorption isotherm can be obtained from Eqs. (43) and (p2[). The result for 6, that is the 
total coverage {6 = c\ + c 2 ), is shown in Fig. 3. We only show 9 as a function of the chemical 
potential (the difference with the representation 6vs. p (oc a) amounts to a logarithmic scale 
change) . 

The staircase form is apparent. The first step corresponds to the filling of the first layer, the 
second one to the filling of the second layer, and so on. It is clear that if we consider more layers 
the structure will persist, but with a progressive collapse of the step size [|], gj. It is of interest to 
determine how the step size for different layers depends on the system parameters, in particular 
on e. The analysis indicates that, for the first layer, A(p, the step's size, as a function of s, scales 
as A(p ~ e u with v ~ 2.0. A similar analysis for the second step indicates that v pa 1.9, showing 
a decrease in step's size, in qualitative agreement with known data [Q, 

It is also possible to analyze the behaviour of the order parameters (that is Aq = cf — c~) 
near the critical point, considering their dependence on «j (that are related with the partial 
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pressure, i.e. cq ~ p). The characterization of such critical behaviour is certainly of interest 
[11]. In our case, from Eqs. (3^,39) and ( ^j49| ) and near the critical point, we have 



with r] = e — 4 (with e ~ 4) ; and 



1 - 1 /2 = I r) i/2 

2 ' 



1/2 _ 11/2 



(so) 



(55) 



(56) 



where rf = e — £q (with eo ~ !)• For both layers, we see that the critical behaviour, as could be 
expected, corresponds to the typical mean field one with a critical exponent 1/2. 



6 Density Profiles 

In this section we show how to get the form of the density profile of a propagating interface 
with constant velocities x>\ for the first layer and i; 2 for the second layer. They can be obtained 
considering the following change of variables in Eqs. (|17|): ci(Ci); with (j = £ — v\t and c 2 (£ 2 ); 
with C2 = £' — V2T. Where ci(Ci) — ► c% when Q\ — ► 00, ci(Ci) — ► C3 when £1 — > —00, and also 
C2(C2) — > c 5 when £ 2 — > 00 c 2 (C2) — ► c 7 when C2 — * -00 

With the indicated change of variables, the partial differential equations (|l^) transform into 
ordinary differential equations 



dc-i , . d , „ . . dc 



vi 



gM) + — (£>(ci)^), (57) 



where 



= 52(c 2 ) + ^( J D(c 2 )^), (58) 
a(,2 «C,2 «C,2 



flf(ci) = ai(l - ci) - ci exp(-eci) (59) 
D(ci) = l-£ Cl (l-ci), (60) 



for the first layer, and 



g(c 2 ) = a 2 {l - c 2 )ci - c 2 exp(-ecic 2 ), (61) 
D(c 2 ) = l-ec 2 (l-c 2 )ci, (62) 

for the second layer. 

For convenience we assume that the concentration jump is localized at £ = 0. It is worth 
noting that with this change of variables the coverage behaves as 

1 /2 

ci(Ci) = c I - > Ci +0, 

and 

ci(Ci) = cj - gxCi , Ci -> -0, 
with g some constant. Similarly, for the second layer 

C2(C2) = c^T - axC 2 /2 , C2 -> +0, 



and 

.1/2 



C2(C2) = cj" - axC 2 , C: 



2 
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with a some constant. 

The mass balance requires that the flux be the same on both sides of the profiles moving 
with velocities v% and v 2 Respectively. Hence, we consider the equation 

D( C1 )^- -D{cC& = -«x(cf-0. (63) 

D(c 2 )^ - D &)^ = -M4-C2)- (64) 

We multiply by D(ci)^, with i = 1,2, respectively. After integrating from — oo to oo; we 
obtain for the velocity of the first layer 

= G{c 3 )-G(ct)-G( Cl ) + G { c+) 

while for the second layer we have 

G(c 7 )-G(4)-G(c 5 ) + G { 4) 

that in general is different from the velocity in the first layer. Note that the propagation velocity 
goes to zero if the numerator of the previous expressions goes to zero. This allows us to represent 
the Maxwell conditions for the coexistence in stationary conditions, for phases subject to first 
order transitions 

G(c 3 )-G( Cl ) = G(c+)-G(c7) (67) 
G(c 7 )-G(c 5 ) = G(c+)-G(c 2 -) (68) 

The profile in the first layer is described by the following integral 

& D(c)dc 



for c x > c > ci , and 



f(c) 



C(c) 



for C3 > c > c\. 

For the second layer we have 

e'(c) 

for > c > C5, and 



[2x(G( Cl ) - G(c))P 
"[2x(G(c 3 )-G(c))]0-5 

[2x(G(c 5 ) - G(e))P 



[2x(G(c 7 ) - G(c))F 

for C7 > c > cj. For C = the coverage jumps at c\ = c\ at ci = Cj , or in the second layer for 
ci + c 2 = c+ + c 2 at ci + c 2 = c\ + cj. 

Figure 5 shows an example of a typical form of the density profiles. It is worth remarking 
here that the form of these profiles is reminiscent of those found in the problem of layering in 
wetting [15]. However, such an analogy requires further study. 
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7 Conclusions 



We have introduced a simple model able to describe qualitatively multilayer adsorption and the 
formation of ordered structures. It is an extension of an early model describing a system of a 
monolayer of adsorbed particles with lateral attractive interactions leading to the formation of 
interfaces separating phases with different local coverages. However, for sake of simplicity, we 
have focused in the case of two layers. Such a case show evidences indicating the occurrence of a 
"cascade" of first order phase transitions in the successive layers. Such a result is in agreement 
with some experimental measurements Q as well as simulations ||. 

Lateral interactions play a key role in the whole picture. Also the possibility of transfer 
interactions among different layers can strongly affect the results, particularly the position of 
the phase transition and the scaling of the step size Aa as a function of e. The study of this 
exchange effect will be done in a forthcoming work. Here we have found that Aa ~ e u with 
v ~ 2.0 for the first step while v ~ 1.92 for the second. We have also found that the critical 
behaviour corresponds to a mean field one, indicating that the analysis of critical properties 
requires a more detailed study. 

The study of the fronts connecting different stationary states [O], for instance the depen- 
dence of the front velocity on system's parameters, as well as the comparison of isotherm with 
simulations and experiments |[L6fl , is under way. Also, the connection of this problem with the 
determination of more realistic adsorption isotherm as well as with cases involving reactions 
[13, 14 1 or layering problems in wetting [15] will be the subject of further work. 
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8 Appendix 



In this appendix we show how to derive the macroscopic master equation starting from a mi- 
croscopic analysis. To reach such a goal we closely follow the procedure of Ref.||], however with 
the adequate changes for the actual situation. The basic step is to consider that the system 
is composed of a given number of particles enclosed in a box and that this number is so large 
that the local coverage, defined as the ratio between occupied sites to the total number of sites, 
should not suffer an appreciable change due to individual adsorption, desorption, or diffusion 
processes. 



8.1 Microscopic approach to the system's evolution equations 

For a box j containing rij particles, the probability per unit time that an adsorption or desorption 
occurs in the first layer is proportional to the number of sites occupied or free in the first and 
those free in the second layer. Hence, for the first layer we have 



w a {n 3 ) = w a (N - n 3 )(l - n'j/N'), (69) 



w d (n j ) = wjn j (l-n' i /N'). (70) 



Where w a = k a p and Wj = k dt o ex P(Ui / ksT) , and N' is the maximum number of sites. 

Here the transition probability per unit time for transitions between neighbor boxes 
(diffusion) is proportional to rij, the particle number in the box j, and the fraction of free sites 
in the box towards which the transition occurs 



wf = wf(l-n J±1 /N) nj . (71) 

Where 



w f = ve k BT Uj < U j± i, (72) 

wf = v Uj> U m . (73) 
Similarly, for the second layer we can write 

w' a (n' j ) = w' a (N' -n^nj/N, (74) 

w'din'j) = wfn'y (75) 

Where w' a = k' a p, «/• = k' d exp^'/fceT), and N' is the maximum number of sites in the second 
layer (N < N'). 

w'^w'+il-n'^/N')^. (76) 



Where 

u'.-u' j±l 



j 



U' < U' ±l , (77) 



wf = v' U'j > U' j±1 . (78) 

For the first layer we assumed N = N max » 1, a number of particles rij and Xj sites 
with j = 1, • • • ,m. The concentration is defined as Cj = -d-. For the second layer we consider 
N' = J2 n jy with N' < N max and N' >> 1. Here we have a number of particles n'j and Xj sites 

nf- 

with j = 1, • • • , m, and the concentration is defined as c'j = jfr. 

We will try to obtain the master equation from the multidimensional probability distribution 
p(ni , ■ ■ ■ , n m , t) . Such a distribution corresponds to the probability of finding n\ , ■ ■ ■ , n m particles 
in the box located at x±, ■ ■ ■ , x m , at time t. Let us see what we obtain from such an analysis for 
the first and second layers. 
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8.2 Equation of motion for the first layer 

The analysis of the time derivative of the probability distribution p(ni, ■ ■ ■ ,n m ,t) for the first 
layer gives 

M{ 2 },t) = ^t^^m-n j + l)pJ-(N-nM{n j },t)) + 
j'=i 

m (N' — n'-) 
+ T, w r N , 3 l(n ] +l)p+-n J p({n :i },t)] + 
3=1 

m 



i 1 



3=1 



" T,W[2- n ^ 1 ^- 1 ] P ({n 3 },t) + 

3=1 

3=1 
m 

+ ElW %1 7 H MHt| (79) 

3=1 



Remembering that 



N' = J2nj, (80) 

3=1 

the local coverage in the second layer is related to 



n'- 

4 = ( §1 ) 



where iV >> 1. These are the available sites for adsorption on the second layer. Also iV >> 1 
and the local coverage for the first layer is Cj = j^. 

Assuming that Cj changes only slightly with an adsorption or desorption process, we can 
expand 

, 1 dp 1 d 2 p , . 

hence, 

D ± ~ v + LA - JP-] + J_[I^ + 1 Q2 P _ d "P 1 f 83 ) 
P > ~ P N [ d Cj dc j±1 l N 2 [ 2 dc] 2 8 2 c J±1 d Cj dc m 1 ' [ j 

Replacing this into Eq. ( |79| ) we get 

3=1 J 3=1 ^ 

m o 1 m o 2 

" E " C 3)^["^] + ^ E-3 (1 - 4)^2 M " 

3=1 3 3=1 3 

m 3 

~ E ^r[( a 3+l C 3+l + Oj'-lCj-l - 2c,-(7 i )(l - C,')p] - 
3=1 J 
m ^ 

~ E^r[( C 3+l+ C 3-l- 2c 3) fT 3 c 3P] + 
3=1 3 
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1 m d 2 

+ 2N 12 {[(^i+icj+i + CT i-i c i-i - ZcjOj-X 1 - cj) - 

j=l 3 

~(cj+i + Cj-i - 2cj)o-jCj]p} - 

i ™ a r a , . a , 
- jfgs^ls^P-^O + s^C-^-i)- 

-S^T^-CjJKo-jCjp) - 
m 3 

" Z ^r[[7i c i( c i+i - c j-i) + (cj - 1)(7j+iCj+i - Ij-iCj-iM - 

j=l 3 

1 m d 2 

~ wl2 75^2 [[^' c i( c i+ 1 - c i-i) + ( X - c i)(7j+iCi+i - 7j-ic,-i)b] - 
ZJV j=i ac j 

As indicated above, the size of an individual box is much smaller than the characteristic 
length scale of the spatial patterns. Therefore, both coverage, Cj and c'j, do not significantly 
change between neighbor boxes, and we can assume that the coverage c(x) is continuous. Hence 
we consider a continuous version of the equation in terms of c(x). These transformation also 
transform the multidimensional distribution p({cj},t) into the functional p(c(x),t). The evolu- 
tion equations for such a functional become 



,9 f , <5 r/ N^ 2 ((TC) <9 2 C. , . , , 

l °J dX m [{l - C) ^ +aC ^ ]P{C{X) ' t) 



5c(x) 

+ Yj dX ~S? [{1 " C ' )iWa(1 ~ C)+ ^ + " " *0]p(c(ar), t) 

+ i/^W ll27C &- 2(1 " C, -^ WC(:C),i)l 
" ^/ <^;C fcM^ [ fcM (1 - <!(;C))l^rC! ' (<!(;C) • ^, 

- ^/^«S5) )V(1 - C),,(C(IM)) - (85> 

Here we have introduced the parameter /j, = where l a is the typical size of each box. This 
value of fj, give us the number of lattice sites per unit area. . 

The coefficients a and 7 in the previous equation represent functions of the coordinate x 
given by 

3 3 (86) 



and 



Here 



2 



W t ~ W i , s 

13 = 3 2 3 ■ (87) 



wj = ve "bt Uj < U j±1 , 
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and 

wj = v Uj> U j±1 . (89) 

Hence, 



and 



a{ x ) = -[1 + e *flTle- l] 5 (90) 
7 (x) = --[l-e ,f fl Tl8 - l ]N, (91) 



with N indicating the sign of 

Considering the limit l Q — > 0, we get 

lim, o ^ (7(^)io) = A (92) 

where 77 is the diffusion constant 

77 = lim lo ^ {vl 2 o ), (93) 

and 

r , m i , r 9U, 1 D dU . 

hm lo ^oh(x)lo] = hm^- — -^-] = "2^f^- (94) 

Introducing the functional differential operator 

" <95> 

considering the evolution equation in the limit Z — > 0, and using the previous equations we 
obtain 

r)r> f 8 D r)TT rfr 

It = -/^)[(l-c')Ml-c)-^ ] + -(--c(l-c)) + t b + 



Dc dU dc ^ J, 7>c 317,, 

+ ( 1 - c )^(^^t)]p- 



+ — / dxl 2 (x)[(l - c')K(l - c) + u> d c] + 77(1 - 2c) 



K B Tdxdx v dc KbT dx 
D f ~ r) 2 ~ 1 3 ~ 1 r?77 

- 7/^ W ^ (4(l)(1 - c))c -5& ( - 4( ^ ) 7^T c(1 - c >fcl p - < % > 

We now consider those terms inversely proportional to fi and independent of the potential 

U, 

F d = ^jdx[A\x)[\(l-2c)^)-A(x)^[A{x){l-c))c)p. (97) 

This expression can be further transformed into 

77 f .«,, . r l. . <9 2 c <9 2 c, - . cL4(x) dc , <9 2 ~ , , 

= -y^^ ( x ) [- ( l-2 c) ^ + c^] + 2^ ) ^- C -^ ) ^[A ( x) ( l-c)c]b 

= - J d X [A\x)[-— 2 ] + c-Ai- - - c)c + y 2 (l - c)c]]p 



A 4 

,d'A 
■ I dx( 



F d =- dx( — ) 2 c(l-c)p, (99) 
u J ox 
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Fd = ~ 1 1 dxdy ^^^^^-^ x ^ 1 - < x ^ x - y ^ ( 10 °) 

Fd = J^J I dxd yM x )My)J^l2Dc(x)(l - c(x))5(x - y)]p. (101) 



Now we consider in Eq. (|9q ) the terms of order 0(fi 1 ), that are proportional to the gradient 
of the potential U. We find that such terms cancel each other, i.e. 

dx[A 2 (x)[-—c— + (1 - c) — (c— ) + — (A 2 (x)c(l - c)—)]p 



2/iksT J dx dx dx dx dx dx 

D f , dU dc . . d . dU. d . . ,017.. , . 

= wr^ (I)[ -& c fc +(1 - c) fe (c & ) -fe (c(1 - c, ^ )b=a (102) 

Using the previous results the evolution equation can be written as 



^ = - [dx-^[w a (l-c')(l-c)-W d (l-Cy+-^^( C (l-c))^)+Dp; 

dt J 5c(x) ksT dx dx dx' 

Iff,, 5 2 



+ tJJ dxdy wm^ - c ' )[ - a(1 - c) + wcmx - y) + 

The last expression corresponds to the functional Fokker-Planck equation for the probability 
distribution functional p(c(x),t). From the theory of stochastic processes, this Fokker-Planck 
equation is equivalent to a stochastic partial differential Jl7| ] (Langevin-like) equation for the 
fluctuating field c(x,t). It has the form 

rlr T) r) f)TT 

- _ M » ( 1_ C)(1 _ C ')_^( 1 _ C ') + __ [C ( 1 _ C )_J + 

+ + A-K(l - C0(1 - c)]'/a(x,() + 4V<0 - <^)]i/i(l,t) + 

+ -4|[(2^(l-c))^/(x,i)]. (104) 



Here f a (x,t),fd(x,t) and f(x,t), are independent noise sources, with intensity one according to 
the Ito interpretation. Those three noise sources correspond to internal noises associated to ad- 
sorption, desorption and diffusion respectively. When studying the equation in the macroscopic 
limit, such noise sources are neglected. 

8.3 Equation of motion for the second layer 

Next, we analyze the time derivative of the probability density p({n'j},t) for the second layer. 
With similar arguments as before we get 

d P ({n'},t) ,„^ n 



w 



+ E<[K + 1 )^i f -^ / (K>'*)] + 

m n' — 1 r/ — 1 

J'=l 
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e + ma - - (i - + 

Eb>i( j+ v — y(K>^)]- (105) 



We should remember that 



m 

AT' = X>i> (106) 



and the coverage of the second layer is related to 

77.- 

4 = ( 107 ) 

where N' >> 1, N' < N max ; and iV 7 are the available sites on the second layer (that is the 
occupied sites in the first layer). 

As in the case of the first layer we approximate 

~± , 1 dp' 1 d 2 p' 

*' 3 



and 



Replacing into Eq. ( |105| ) we obtain 

MWjht) A , a d [n i - , a ^ rn , 

- e^kvi+^e^^V)- 

5=1 J 5=1 J 

m pi 

- e^[^'( c ; + i- 2 4+4-i)]+ 

5=1 3 
1 m 

+ ^E^[(^ +1 4 +1 + ^ 1 4_ 1 -2cX)(i-4)- 

- ^ e ^[^(i - + - ^) - 2^d - - 

- E ^7^4(4+1 - c i~i) + (4 - 1 )(^'+l4+i - 7j-i4_i)b' - 

j=l UC j 

I m d 2 

~ ^ [7i c j( c i+i " c i-i) + C 1 " c i)(^+i c i+i " lj-l c j-i)]p' ~ 

i=i ac i 

-i 772 q o 
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As for the first layer we consider a uniform coverage implying 

c( Xj )^c(x), (111) 

c'(xj) - c'(x). (112) 

Also 

p({cj},t)^p(c(x),t), (113) 

V'{{^t)^p{d{x),t). (114) 

The continuous version of the equation of motion, assuming the uniform coverage of the 
system as before, is again a functional equation for p(c'(x),t). This equation is 



d P (c'(x),t) 
dt 



+ 27/ ^^[[(l-^c + ^c' + Z^l-c^ 2 ^ 



2/i'7 Sd(x) 2LlK ' oy ' dx 2 

-a'l 2 c^]p({c'(x)},t)} + 

+ 27/ d Wf 7 ' C '£ " " C ^H{c'(,)},t)] - 

- ^/^(iWb 1 ^^ 1 - 40 ^ 1 '^ (115) 

Again as in the first layer case, we have introduced here the parameter ji' = j-, where l a is 
the size of each box. The value of // gives the number of sites per unit of area. 

The coefficients a' y 7' in the last equation differ from those in the first layer by a factor of 
H 1 . Hence, 

a'{ x ) = -[1 + e JC B T ' 8 * 1 ], (116) 



2 
"2" 



7 '( x ) = [1 _ e jr fl TlftJ] ^ (117) 



where N, as before, indicate the sign of 
In the limit l a — > 0, we get 

/im /o _ (7'W^) = (118) 
where D' is the diffusion constant on the second layer, 

D' = lim lo ^{v'l 2 ), (119) 

and 

We will adopt £)' ~ D, as they differ by a constant. 
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We introduce now the differential operator 

Introducing it into the evolution equation and considering the limit l D — ► 0, we get 
I = -/^)[K a (l- C 0-^] + ^(^fc'(l-cO) + 
+ D^]p(c'(x),t) + 

+ T~> J dx B 2 (x)llcw' a (l ~ c) + w ,d c'} + D(l - 2c')^ - 
Dd dUdd . d . Dd dU^ , . . 

+ ( 1 - c )^7(^7^r)Mc'(x),0- 



dx v ' dd K K B T dx 

D [ , , d 2 



^ 7 y ^[5(^^(5(^(1 -c')) 

1 r) ~ 1 dTT 

-2& (BW, »^ (1 -% M</|lM) ' <122) 

As for the previous first layer case, we consider the terms inversely proportional to (j! and 
those independent of the potential U. 



This expression can be transformed into 



D r 1 d 2 d ~ f) 2 „ 

= -y^[B 2 (x)[-(l-2 C 0^]-i?(x)^[i?(x)(l-c')c']Kc'(^),i)- (123) 



w # /" , rA 2/ s r 1 / A aV .dV, ,dB(x)dd , 

-B(x)^[B(x)(l-d)d]]p 

D f , rfi 2 / Nri^c'n ,dB 2 (x)dd 1 d 2 rft2 , Wl 
= 7 7 ^ ix)[ 2d^ ] + C ^ " 2 dx~ 2[B W(1 " c )c 

+ [ H ]2(1 " C ' )C ' M 
D f S2 1 5 2 c' ldB{x)dd dB 



7 ^Wbd + o-^*r + hjrK 1 " C ')c>, (124) 



\n j 2 dx 2 2 dx dx dx 

F 'd = ^j dx{ ^) 2 ^ 1 ~ c >' ( 125 ) 

F ' d = ~< 1 1 dxdy( ^^ )( ^^ )c ' (x)(1 " C ' (x))5(x " (126) 

F ' d= 2j'S J dxd vB^)B(y)^[2Dd(x)(l - c'(x))5(x - y)]p. (127) 



Now we consider in Eq. (|122| ) the terms of order 0(/i x ), that are proportional to the 
gradient of the potential U. We find that such terms cancel each other, i.e. 

= w/ <ix|B(x)| -^ +(1 -^ ) & ( ^ )l 
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2/i'fcsT J dx dx dx dx 

Taking into account our previous results, the evolution equation becomes 



dp _ _ f , 8 
dt J 8c' (x) 



f dx—^—[w' a (l - c')c - w' d c' 
J oc[x) 

D d , , S .8U, n 5V. 



+ v J J dxdy ^wm m ~ c >'* c + w ' d ^ x ~ y ^ 

+ a5^ (2jDc/(i " c ' ),5(x " 2/))b - (i29) 

The last equation is the functional Fokker-Planck equation for the probability distribution 
p(c'(x),t) for the second layer. Similarly to the discussion for the first layer, we can obtain 
the related stochastic partial differential equation 

dc' ,„ /N ,a , D d . ,,911, 



,#V 1 „ ,i i 

dx 2 



+ + ^[w' a (l - c')c]$f' a (x,t) + ^[«/ d </]3/dOM) + 

[12 /J, 2 



-l-A[(2D c '(l- c '))t/'(x,t)]. (130) 



Here, as before, f^(x,t),f' d (x,t) and f'(x,t), are independent noise sources, with intensity one 
according to the Ito interpretation. These three noise sources correspond to internal noises 
associated to adsorption, desorption and diffusion respectively, but now in the second layer. 
Also, when studying the equation in the macroscopic limit, such noise sources are neglected. 
The previous analysis leads us to Eqs.(17) and (p^). 
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Figure 1: Sketch of the different processes included in the model. 



Figure 2: Phase diagrams in the planes e vs. at. a) first layer, b) second layer. Region I 
corresponds to the dilute (vapor-like) phase, region III to the dense (liquid-like) phase, while II 
corresponds to the region of phase coexistence. 



Figure 3: Typical form of the isotherm. We depict the total coverage 6(= c\ + c-i) vs. tp, the 
chemical potential. 



Figure 4: Here we depict the relation between the step size and the potential parameter e: Aip 
vs. e, yielding a linear dependence in logarithmic scale. The result here corresponds to the 1st 
step. 



Figure 5: Typical form of the density profile. Here we have used e = 5.6, ot\ = 0.08 and 
a 2 = 0.096 
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